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The  goal  of  this  research  project  is  to  measure  local  grain  size  and 
concentration  deep  in  the  interior  of  a  medium  ultrasonically.  It  will  be 
shown  that  information  on  the  scatterers*  statistics  can  be  obtained  if 
either  the  field  structure  or  the  scatterer  density  exhibits  gradients. 

It  will  be  shown  theoretically/  and  by  numerical  computation/  that  the 
received  echo  from  a  configuration  of  piston  transmitter- and  a  point 
receiver  in  the  center  of  the  piston  contains  a  coherent  component  which 
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..should  a UQ'i  concentration  estimation.  This  report  will  also  show 
theoretically  and  confirm  experimentally  that  a  gradient  in  scattering 
concentration  will  return  a  coherent  echo  from  whose  degree  of  coherence 
th e  scatterer  concentration  can  be  estimated. 

Ue  expect  to  apply  these  results  soon  to  the  estimation  of  graijr^ 
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A.  INTRODUCTION 


I.  Project  Objectives 

The  goal  of  the  research  described  in  this  report  was  to  find  techniques 

for  estimating  local  values  of  the  grain  density  in  metallic  and  ceramic  media, 
using  non-destructive  ultrasonic  examination.  Although  the  goal  of  non-destruc¬ 
tive  ultrasonic  local  grain  size  estimation  deep  in  the  interior  of  a  medium 
has  been  a  subject  of  widespread  interest  for  a  number  of  years,  no  technique 
for  accomplishing  this  had  been  found  at  the  outset  of  this  project.  The 
techniques  for  the  ultrasonic  determination  of  grain  size,  which  are  based  on 
the  variation  of  ultrasonic  attenuation  with  grain  size,  do  not  permit  the 
determination  of  grain  size  at  a  point  inside  a  test  sample,  but  merely  allow 
the  determination  of  the  grain  density  averaged  over  the  length  of  the  sound 
beam  from  the  sample  surface  to  the  interior  point  of  measurement. 

II.  Research  Achievements 

During  the  period  covered  by  this  contract,  it  has  been  possible  to  develop 
a  technique  for  analyzing  backscattered  echoes  which  has  been  shown  to  allow 
the  determination  of  grain  size  at  medium  boundaries  and  which  shows  great 
promise  of  permitting  grain  size  determination  at  interior  points  also. 

The  technique  for  surface  grain  size  determination  is  based  on  a  theory 
described  in  Section  B-II  which  shows  that  the  backscattered  echo  from  an 
ensemble  of  grain-like  scatterers  contains  both  coherent  and  incoherent  echoes, 
from  whose  ratio  the  surface  grain  density  can  be  estimated.  Section  B-III 
describes  experiments  on  sponges  which  confirm  the  results  of  the  theory,  and 
which  suggest  that  the  method  should  be  applicable  to  metals  and  ceramics. 

Section  C-I  extends  the  theory  and  shows  that  for  focussed  ultrasound 
fields  of  sufficiently  low  symmetry,  illuminating  grainy  media,  it  should 
again  be  possible  to  obtain  coherent  and  incoherent  echoes,  from  whose  ratio 
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the  grain  size  can  be  determined.  The  section  contains  a  calculation  of  the 
coherent  echo  for  a  specific  transducer  geometry  which  demonstrates  the  effect, 
and  whose  program  is  reproduced  in  the  appendix  (section  F). 

Metallurgical  preparation  of  test  specimens  for  this  research  is  described 
in  section  D. 


III.  Conclusions 

This  research  has  introduced  and  provided  the  theoretical  foundations  of 
a  novel  technique  for  non-destructively  estimating  the  density  of  scatterer 
ensembles  from  the  ratio  of  their  coherent  to  their  incoherent  ultrasonic 
echoes.  It  was  shown  that  the  coherent  echoes  required  for  this  process  are 
produced  either  by  gradients  in  scatterer  density,  or  by  field  foci  with  suf¬ 
ficiently  low  symmetry.  The  effectiveness  of  the  technique  for  estimating 
surface  scatterer  density  was  demonstrated  experimentally  on  sponges.  In 
future  research  we  will  attempt  to  demonstrate  that  the  technique  using  focussed 
beams  can  be  used  to  estimate  the  grain  density  in  the  interior  of  metals 
and  ceramics  of  practical  importance. 


I.  Introduction 

It  Is  known  that  a  uniformly  illuminated  infinite  ensemble  of  random  scat- 
terers  produces  an  incoherent  echo,  (defined  as  having  zero  mean  at  any  given 
range  delay),  whose  power  is  proportional  to  Np2.  Here  N  is  the  scatterer 
volume  density,  p  the  scattering  coefficient,  and  represents  ensemble  aver¬ 
age.  It  has  been  predicted,  although  apparently  never  experimentally  verified, 
that  a  coherent  echo  should  be  generated  by  random  scatterers  either  when  these 
exhibit  a  density  gradient  *■,  or  when  they  occupy  a  bounded  region2 *3.  We  have 
now  confirmed  these  results  by  experiment  and  show  that  an  estimate  of  the 
density  can  be  extracted  from  the  measured  ratio  of  the  coherent  and  incoherent 
echoes . 


II .  Theory 

Let  the  echo  detected  by  a  sonic  or  electromagnetic  receiver  at  time  t 
after  the  transmission  of  a  burst,  due  to  a  single  point  scatterer  at  r  be 
written  as 


Er(r.t)  ■  PG(r,t)  (1) 

where  G(r,t)  is  defined  as  the  system  impulse  response.  When  many  scatterers 
are  illuminated,  the  echo  at  time  t  after  transmission  can  be  written 

E(t/n)  -  l  P1G(r1,t)  (2) 

i-1 

where  n  is  the  number  of  scatterers  in  the  "range  cell",  defined  as  that  region 
of  space  from  which  echoes  are  received  at  time  t  after  the  start  of  transmis¬ 
sion.  We  assume  that  this  region  has  volume  V,  and  that  it  contains  n  randomly 
positioned  scatterers  at  locations  rj,  r2,  of  strengths  pi,  P2«».  pn» 

The  instantaneous  power  at  time  t  is 


P(tjn)  -  l  l  pj  G(r1,t)  G  (rj,t) 
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This  can  be  separated  into  two  parts 

P(t|n)  -  l  pf  |G(r±,t) |2  +  I  I  Pi  P4  G(r^,t)  G*(r1,t)  (4) 

i«l  1  j  J  J 

From  eqt .  (4)  the  average  power  at  time  t  after  transmission,  given  that 
there  are  exactly  n  scatterers  in  the  range  cell,  can  be  written 

P( 1 1 n)  «  n  p£|G(ri,t)|2  +  n(n-l)  Pj  Gfr^t)  G(rj,t)  (5) 

where  n  is  the  number  of  scatterers  in  range-cell  volume  V.  Note  that  the  aver¬ 
age  P(t | n)  is  on  the  random  variable  r^  and  P? 1 1 n)  is  the  average  over  all  pos¬ 
sible  configurations  of  the  r^’s  in  the  range  cell*  Assuming  that  the  scatter- 
ering  coefficients  Pi  are  uncorrelated  with  each  other  and  that  ri  is  uncor- 
related  with  the  rj  for  i*j,  we  can  rewrite  eqt.  (5)  as 

P(t|n)  ■  n  ~pZ  f  |G(r,t)|2  p(r)dv  +  n(n-l)  p2  |  /  G(r,t)p(r)dv j 2  (6) 

V  V 

where  p(r)  is  the  probability  of  finding  one  individual  scatterer  in  the  vol- 
ume  element  at  r,  and  where  the  integrals  are  taken  over  the  volume  V  of  the 
range  cell.  We  assume  the  scatterers  to  be  uniformly  distributed  in  the  range 
cell.  So  we  may  write 

P(r)  -  I 
V 

Using  these  two  relations  we  can  write  the  ensemble  averaged  power  P(t |n) 

P( 1 1 n)  -Bp2/  |G(r,t)  |2dv  +  *>(n-l)  |  /  G(r,t)dv|2  (7) 

v  V  V2  V 

To  find  the  average  power  p( t)  we  have  to  evaluate 


pTtT  "  l  P(t|n)  p(n) 


(8) 


from  7  and  8 


p(t)  -  /  |G(r,t) | 2dv  l  np(n)  +  11  | /  G(r,t)dv|2  J  n(n-l)p(n)  (9) 

v  V  n-1  V2  V  n-1 

where  n  is  Poisson  distributed  with  n  as  its  average.  Replacing 

00 

I  np(n)  -  n 
n-1 

and  for  Poisson 

00 

J  n(n-l)  p(n)  -  n2 
n-1 

we  get 

pOO  -  N  p*  /  |Gr,t)  |  2dv  +  N2  p2  |  /  G(£,t)dv|2  (10) 

V  V 

where  N  -  £ 

V 

We  see  that  the  first  term  on  the  r.h.s  of  eqt.  (10)  is  proportional  to  Np7, 
i.e.,  to  the  average  of  the  sum  of  the  individual  backscattered  powers.  We  will 
therefore  refer  to  it  as  the  incoherent  power  Pine*  The  second  term  is  propor¬ 
tional  to  (Np)2,  i.e.  to  the  square  of  the  sum  of  the  backscattered  amplitudes. 

We  will  therefore  describe  it  as  the  coherent  backscattered  power,  Pcoh* 
We  now  show  that  Pcoh  equals  the  magnitude  squared  of  the  averaged  echo, 
and  that  P^nc  equals  the  variance  of  this  echo. 

From  eqt.  (2) , 


E(t)  -  l  E(t | n)  p(n)  -  Np  /  G(r,t)dv 
n-1  V 


(ID 


So 


Pcoh  "  iEl2  “  n2  P2  1/  G<f.Odv|2 


(12) 


2 
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The  fact  that  Pinc  is  the  variance  of  E(t)  follows  immediately,  since 
Thus  from  eqts.  (7)  and  (9), 


-  IEl7  -  | E | 2 


°E  "  pinc  "  "P2  /  !G(r,t)|2dv  (13) 

V 

Let  us  now  examine  the  conditions  under  which  the  coherent  backscattered 

echo  can  exist.  From  eqt.  (7)  the  ratio  between  the  coherent  and  incoherent 

echo  powers  is  pcoh/pinc  “  NF»  where  F  -  p2  |  /  G(r , t)dv 1 1/'pz  /|G(r,t) | 2dv. 

For  p  /p  of  order  unity,  the  coherent  echo  will  only  be  appreciable  if 

both  the  scatterer  density  N  and  the  function  F  are  sufficiently  large.  To 

examine  the  magnitude  of  this  function  we  start  from  the  observation  that  from 

00 

-*•  ;  * 

the  definition  of  G(r,t)  given  in  eqt.  (1),  it  must  follow  that  J  G(r,t)dt  =  0 

o 

since  an  ultrasound  echo  can  not  contain  a  d.c.  component.  We  then  see  that 
although  the  incoherent  echo  power  which  is  proportional  to  f  |G(r,t)|2dv  will 
always  be  greater  than  zero,  the  same  is  not  necessarily  true  for  the  coherent 
power  term  which  is  proportional  to  |G(r,t)dv|2.  For  instance  for  a  plane  wave 
travelling  along  the  z  axis,  the  field  function  can  be  written 

G(r,t)  -  I(t-lz)  ei<ut“2ltz) 


Thus  from  eqt.  (9) 

E(t)  ■  pN  fdy  f dx  /l( t-2^)e^(u,t_2^z^dz 

c 

*  > 

which  is  equal  to  zero  since  /  G(r,t)dt  *  0 

o 

Likewise  for  a  point-transmitter-receiver 

+  ♦  ♦ 
G(r,t)  -  [I(t-2r/c)/r2]e1(“t_2k,r) 


E(t)  =»  pN  f  [I(t-2r/c)/r2]e^(u,t-2k*r^  •  4irr2dr 


which  again  reduces  to  zero.  Thus  we  see  that  both  for  the  case  of  plane  waves 

and  for  the  case  of  waves  emitted  from  a  point  transmitter/receiver ,  the  beam 

-► 

geometry  has  such  a  high  degree  of  symmetry,  that  the  volume  integral  f  G(r,t)dv 

> 

is  found  to  be  equal  to  f  G(r,t)dt,  which  has  to  be  zero.  Thus  no  coherent 


signal  exists  in  these  two  geometries  for  uniformly  distributed  scatterers. 


Density  Gradient 

_  *► 

If  the  scatterers  exhibit  a  certain  density  profile  N(r),  the  algebra  in¬ 
troduced  earlier  becomes  somewhat  more  involved,  setting  p(r)  to  be  •  A  the 

n  v 

average  power  at  a  certain  range  delay  t  can  be  shown  to  be, 

p(t7  «  f  N(r)|G(r,t)|2dv  +  p2  | /  N(r)  G(r,t)dv|2  (14) 

V  V 


with 

E(t)  -  p  /  N(r)G(r,t)dv 
V 


we  can  now  investigate  the  properties  of  the  returned  echo  from  the  boundary 
between  2  regions  with  different  scattering  densities.  For  simplicity  we  as- 
sume  the  imaginary  boundary  to  be  planer.  Replacing  r  by  z  we  can  write, 


N(z)  -  Nj_ 
-  N2 


z  >  zo 

Z  <  Z0 
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(15) 


For  a  plane  wave  G(z,t)  <*  b(  t-AiDexpl-ioKt-Zz/C]  ;  and  the  returned  echo  from 

C 

the  boundary  becomes 

o  °> 

E(t0)  *  Ni  pY  /  b(t)exp[-(Dt]dt  +  N2  P2  /  b( t)exp[-u>t)dt  (16) 

—oo  O 

PI  and  p2  are  the  reflection  coefficients  of  the  scatterers  in  region  1  and  2 
respectively.  Note  that  the  volume  integral  In  14  becomes  a  time  integral  by 
change  of  variables.  If  the  range  cell  is  situated  such  that  the  boundary  is 
closely  in  its  center,  the  variance  of  E(t0)  can  be  shown  to  be, 

00 

o2  s  A[0^  +  o|]  *  [NiP^  +  /  |l»(t)|2dt  (17) 

2  o 

2  2 

where  and  a ^  are  the  variances  of  the  returned  echo  from  region  1  and  region 


2  respectively. 


III.  Experimental  Results 

According  to  the  results  obtained  earlier  we  should  be  able  to  observe  a 
coherent  effect  from  an  echo  returned  from  the  boundary  between  2  different 
scattering  density  regions. 

One  of  the  easiest  ways  to  obtain  scatterers  immersed  in  water  and  have  a 
sharp  boundary  is  to  use  sponges.  The  sponge  can  be  cut  to  produce  a  sharp 
boundary.  The  sponge  is  then  Immersed  in  water  and  the  complex  structure  of 
the  sponge  fibers  can  be  considered  as  randomly  distributed  scatterers  immersed 
in  water.  Furthermore,  we  can  clamp  together  two  sponges  with  different  scat¬ 
tering  properties  and  compose  a  sharp  boundary  between  2  different  media,  or  we 
can  use  only  one  sponge  for  which  we  have  2  media  with  Ni*0  In  region  1  and  N2 
in  region  2.  Figs,  la-ld  show  pictures  of  sponges  used  in  the  experiments. 
One  can  see  that  sponge  A  has  the  finest  grain  structure,  sponge  B  has  a  larger 
honeycomb  structure,  sponge  C  is  more  dilute,  and  sponge  D  is  very  dilute  in 
comparison  to  sponges  A,  B  and  C.  In  fig.  2  we  see  sponges  B  and  C  side  by 
side.  One  can  see  that  the  boundary  is  sharp  in  comparison  to  a  wavelength 
(frequency  of  2.25  MHz).  One  can  also  see  again  that  the  honeycomb  structure 
of  sponge  B  is  finer  than  that  of  sponge  C.  In  fig.  3  we  see  the  echo  returned 
from  sponge  A.  One  can  see  that  the  echo  from  the  boundary  is  highly  coherent, 

i.e.,  the  echoes  do  not  change  phase  and  are  not  as  random  in  comparison  to  the 

echo  from  Inside  the  sponge.  These  echoes  are  incoherent  in  the  sense  that  for 
a  certain  range  delay  the  amplitude  is  random  with  zero  mean.*  From  the  theory 
introduced  earlier  we  know  that  the  ratio  between  the  average  echo  squared  from 

*Note  that  the  echoes  from  close  vicinity  to  the  first  echo  are  somewhat  weaker 
than  the  echoes  from  the  echoes  deep  inside  the  sponge  (eventually  the  echoes 
are  weakening  due  to  attenuation).  The  reason  stems  from  the  fact  that  the 
transducer  surface  is  not  perfectly  parallel  to  the  sponge  surface  and  thus 

it  takes  more  time  for  the  Incoherent  term  to  fully  develop.  The  slight  angle 

change  has  practically  no  effect  on  the  coherent  term  magnitude.  So  the 
estimation  of  the  incoherent  term  should  be  done  further  away  from  the  first 
echo. 


the  boundary  and  the  average  power  from  inside  the  sponge,  neglecting  attenua¬ 
tion,  is  proportional  to  the  scattering  density.  So,  practically,  one  can 
estimate  the  sponge  density  from  the  results  obtained.  Ve  will  not  attempt  here 
to  estimate  the  sponge  density,  but  will  show  qualitatively  that  the  experi¬ 
mental  results  behave  according  to  the  theory.  Figs.  4A,B,C,  and  D  represent 
the  echoes  received  from  sponges  A,  B,  C,  and  D  respectively.  Each  of  the 
pictures  is  composed  of  four  pictures,  one  on  top  of  the  other,  from  four 
different  locations  in  the  sponge,  with  identical  distances  from  the  transducer 
to  the  sponge.  Fig.  4A  is  similar  to  fig*  3,  with  the  traces  one  on  top  of  the 
other.  One  can  see  that  the  coherent  component  of  the  echo  from  the  boundary 
is  relatively  large  in  comparison  to  the  echo  returned  from  within  the  sponge. 
Sponge  D,  with  the  largest  honeycomb  structure,  has  practically  zero  coherent 
effect.  The  honeycomb  structure  of  sponges  B  and  C  is  similar  in  shape,  so  we 
can  asstime  that  p  /p7  is  also  similar.  The  integrals  in  eqt.  10  are  also 
similar  for  the  two  sponges,  as  they  involve  only  transducer  parameters,  which 
were  the  same  for  the  experiments.  So  the  ratio  between  the  coherent  power 
term  and  the  incoherent  term  should  be  proportional  to  the  scattering  density. 
Evaluating  the  coherent  term  from  the  average  of  the  first  echo  from  the  bound¬ 
ary  and  the  power  from  inside  the  sponge  (the  incoherent  term),  we  get  NF  to  be 
about  2.1  from  sponge  B  and  about  0.8  from  sponge  C,  which  implies  a  density 
ratio  of  about  2.6.  The  actual  scattering  density  ratio  that  can  roughly  be 
established  from  the  micrographs  is  about  2.2,  which  is  in  close  agreement,  con¬ 
sidering  the  fact  that  we  used  only  four  sample  points  and  no  special  arrange¬ 
ments  for  producing  extremely  smooth  surface  and  keeping  the  sponge  surface 
parallel  to  the  transducer  surface.  Sponge  D  has  much  lower  density,  too  low 
to  establish  a  coherent  term  out  of  only  4  sample  points.  In  fig.  S  the  echo 
from  sponge  D  clamped  to  sponge  B  is  shown.  The  coherent  effect  is  clearly 
seen  along  the  center  vertical  line  of  the  picture.  We  investigated  also  the 


angle  dependence  of  Che  reflected  echo*  The  sponge  under  test  was  chosen 
Co  be  sponge  B.  One  can  see  In  fig*  6  that  Che  coherent  effect  is  highly 
angle-sensitive,  while  the  echo  reflected  from  within  the  sponge  is  not  angle- 
dependent*  The  behavior  of  the  coherent  component  is  very  much  like  specular 
reflection  (not  proved  in  this  report);  thus  its  behavior  resembles  specular 
reflection.  The  reflection  from  within  the  sponge  is  independent  of  the  angle 
as  the  echo  is  independent  of  the  boundary  region. 

Conclusions 

He  have  pointed  out  earlier  that  since  neither  an  electromagnetic  nor  an 
acoustic  transducer  can  transmit  d*c.  signals,  the  time  integral  of  all  trans¬ 
mitted  electromagnetic  or  acoustic  signal  amplitudes  is  zero*  This  implies 
that  the  echoes  of  such  signals  from  constant  density  random  scatterers  are 
purely  incoherent  for  either  a  plane  wave  or  a  point  source  emitter/receiver  * 
In  the  presence  of  strong  density  gradients,  however,  such  as  those  which  occur 
near  a  boundary,  we  have  demonstrated  the  existence  of  coherent  echo  components, 
defined  as  echoes  with  a  finite  ensemble  average  at  a  specific  range  delay. 
Measuring  the  ratio  between  the  coherent  echo  and  the  echo  variance  should  per- 
mit  estimation  of  the  quantity  Np  / p  ,  which  in  the  case  of  metals,  for  example, 
is  related  to  the  grain  size. 


K>  n  \  v*  ~  /  K'f  K  *  - 

v>i  &  ■ 

't  -tXM? 

hS2.  V* 


h  *  •  i  '  A  1 

r-\  **- 


Sponges  B  and  C  side  by  side.  It  is  seen  that  the  honeycomb  structure 
of  sponge  B  Is  finer  than  that  of  sponge  C. 
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Four  different  echoes  from  sponge  B.  Note  the  coherent  echo  component 
from  the  sponge  boundary  and  the  Incoherent  component  from  within  the 
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Fig.  4a.  Four  traces  (one  on  top  of  the  other)  of  the  reflected  echo  from 
sponge  A.  Note  the  large  coherent  component  of  the  echo  from  the 

boundary. 

Fig.  4b.  Sponge  B.  Note  that  the  coherent  component  is  not  as  large  in  com¬ 
parison  to  the  echo  from  within  the  sponge  as  for  fig.  4A. 

Fig.  4c.  Sponge  C.  Note  that  the  coherent  component  is  smaller  in  comparison 

to  the  echo  from  within  the  sponge  than  in  fig.  4B. 

Fig.  4d .  Sponge  D.  The  coherent  component  was  found  to  be  negligible  in 

repetitive  experiments. 


The  received  echo  from  a  2  layer  sponge  complex  composed  of  sponges  C 
and  B.  The  coherent  component  can  be  seen  on  the  central  vertical 
line.  Also  note  that  the  power  reflected  from  sponge  C  Is  less  than 
that  reflected  from  sponge  B. 


Angle  dependence  of  the  coherent  component  (sponge  B) .  Note  that 
while  the  coherent  term  Is  highly  sensitive  to  angle  variations, 
the  echo  from  within  the  sponge  is  angle  Insensitive. 


I.  Introduction 


The  technique  introduced  in  the  last  chapter  enables  us  to  estimate  the 
scattering  density  near  a  boundary,  or  when  a  sharp  scattering  density  exists. 
However,  if  we  need  to  estimate  uniform  scattering  density  inside  the  metal, 
this  technique  would  clearly  fail.  However,  if  we  introduce  a  sharp  field 
gradient  in  the  range  of  interest,  a  coherent  echo  will  still  be  generated,  as 
can  be  learned  from  eqt.  (10)  of  the  preceding  chapter.  In  this  section  we 
investigate  the  coherent  echo  created  by  a  field  gradient  expected  from  a 
simple  geometry  of  transmitter/receiver.  The  transmitter  is  a  circular  disk, 
and  the  receiver  is  a  point  transducer  at  its  center.  We  show  that  the  configu¬ 
ration  leads  to  a  coherent  echo  from  a  half-space  filled  with  inhomogeneities. 

II.  Field  Requirements  for  Coherence 

Let  us  now  examine  the  conditions  under  which  the  coherent  backscattered 
echo  can  exist.  From  eqt.  (7)  the  ratio  between  the  coherent  and  incoherent 
echo  powers  is  *  NF  p^/p2,  where  F  m 

/|G(f,t)|2dv 

Since  p^/p2  is  usually  of  order  unity,  the  coherent  echo  will  be  appreci¬ 
able  only  if  both  the  scatterer  density  N  and  the  function  F  are  sufficiently 
large.  To  examine  the  magnitude  of  this  function,  we  start  from  the  observation 
that  from  the  definition  of  G(£,t)  given  in  eqt.  (1),  it  must  follow  that  /“ 
G(f,t)dt  *  0,  since  an  ultrasound  echo  can  not  contain  a  d.c.  component.  We  see 
therefore  that  although  the  incoherent  echo  power  which  is  proportional  to 
||G(f,t)|2  p(£)dv  will  always  be  greater  than  zero,  the  same  Is  not  neces¬ 
sarily  true  for  the  coherent  power  term  which  is  proportional  to  | /G(f,t)p(f)dv| 
For  Instance  for  a  plane  wave  travelling  along  the  z  axis,  the  field  function 
can  be  written 

G(f,t)  -  I(t-2z/c)ei(<*»t-2kz) 


Thus  from  eqt.  (8) 
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E(t)  -  pN  jdy  fix  J I(t-2*)e1(alt"2kz>dz 

oo 

which  is  equal  to  zero  since  /G(f,t)dt  -  0. 

o 

Likewise  for  a  point  transmitter-receiver 
G(?,t)  -  [I(t-2r/c)/r2]e1<“t“2k,r> 

and 

♦  ■* 

E(t7  -  pN  /  [I(t-2r/c)/r2le1<ut_2k*r)  •  4Trr2dr 
o 

which  again  reduces  to  zero.  Thus  we  see  that  both  for  the  case  of  plane  waves 

and  for  the  case  of  waves  emitted  from  a  point  transmitter/receiver ,  the  beam 

+■ 

geometry  has  such  a  high  degree  of  symmetry  that  the  volume  Integral  f  G(r,t)dv 
is  found  to  be  equal  to  jG(f,t)*t  which  has  to  be  zero.  Thus  no  coherent  signal 
exists  in  these  two  geometries  for  uniformly  distributed  scatterers. 

The  above  calculations  suggest  that  a  non-zero  coherent  echo  would  be 
produced  with  separate  transmitters  and  receivers,  preferably  having  opposite 
magnitudes  of  field  gradient  in  the  region  of  maximum  intensity.  Such  a  geom¬ 
etry  would  be  provided,  for  example,  by  a  focussed  transmitter/point  receiver 
pair.  For  reasons  of  computational  simplicity  we  have  chosen  to  examine  the 
case  of  an  unfocussed  flat  circular  piston  transmitter  set  in  a  rigid  baffle, 
with  a  point  receiver  at  its  center,  (see  Fig.  7).  For  this  type  of  transmitter 
Arditi  et  al*  have  shown  that  for  a  velocity  normal  to  the  transducer  face  v(t), 
the  transmitted  sound  pressure  at  point  r_  can  be  written 

p(r,t)  -  -  z  3v/3t  *  ht(r,t)  (11) 

where  e  is  the  density  of  the  medium  and  the  transmitter  impulse  response  is 

ht(r,t)  -  1/2 w  /  6(t  -  r'/c)/r'  dS 
S 
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which  has  been  calculated  fcr  the  flat  disk  piston  transducer  In  a  rigid  baffle 
by  Oberhettinger^. 

The  echo  due  to  a  single  scatterer  at  t_  for  which  p  is  unity,  can  then  be 
written 

G(r,t)  »  p(r,t)  *  hr(r,t)  (12) 

♦ 

where  hr(r,t)  is  the  receiver  Impulse  response  which  for  a  point  receiver  at 
the  origin  becomes 

hr(r,t)  -  5(t  -  r/c)/2irr  (13) 

Combining  Eqts.  (11),  (12)  and  (13)  gives  for  the  echo  of  a  single  scatterer 
at  f, 

G(r,t)  *  -  e/2n  3v/3t  *  [1/r  ht(r,  t  -  r/c)]  (14) 

+  ♦ 

which  reduces  to  G5(r,t)  -  -  e/(2irr)ht(r ,  t  -  r/c)  if  3v/3t  is  taken  to  be  an 

impulse  at  time  zero. 

-*■ 

Integrating  G$(r,t)  over  space  provides  us,  as  shown  by  Eq.  (8)  with 
the  spatially  averaged  coherent  echo  E$(t)  resulting  from  3v/ 3t  ■  6(t)  for 
uniform  average  density  scatterers  for  which  pn  is  unity.  This  integral  is 
plotted  in  Fig.  8  for  t  >  a/c.  It  can  be  seen  that  as  t  approaches  infinity 
Eg(t)  approaches  the  constant  a  c/4ir,  which  corresponds  to  the  echo  that  would 
be  obtained  from  a  point  source/receiver  of  strength  a^,  illuminating  a  half 
space.  By  analysing  the  case  of  a  point  source  in  an  infinite  plane,  it  can 
be  shown  that  near  the  origin  E g(t)  +  c^/4  t^  corresponding  to  the  analytically 
obtainable  solution  for  an  infinite  area  transmitting  transducer. 

For  a  transducer  acceleration  v(t)  of  finite  length,  the  spatially  averaged 
coherent  echo  will  be  given  by 


E(t)  -  v(t)  *  Ej( t)  . 


(15) 


This  quantity  has  been  calculated  for  a  transducer  surface  displacement 
consisting  of  m  sinusoids  i.e.  for  u  *  sin  u>t  for  0  <  t  <  2*m/u),  where  m  is  an 
Integer,  and  with  the  simiplifying  assumption  that  we  can  represent  the  spatial¬ 
ly  averaged  coherent  echo  Eg(t)  resulting  from  impulse  acceleration  and  shown 
in  Fig.  8,  as 

E »(t)  -  a2[K  e'at/a  +  £_]  for  t  >  a/c  (16) 

4* 

with  K  »  30,000  cm/sec  and  a/a  ■  5000Hz. 

Substituting  eq.  (16)  into  eq.  (15)  gives,  after  some  higher  algebra, 

2mira 

1(0  -  a  K  m3  [1  -  e  au)  ]  e-i*/*)1  for  t  >  i  >  . 

(a/a)^+m2  c  “  (17) 

We  see  that  for  the  geometry  analyzed  here,  of  a  disk/point-source  transducer 
pair,  the  maximum  coherent  echo  occurs  immediately  after  the  end  of  the  trans¬ 
mitted  signal  and  therefore  comes  from  a  region  close  to  the  point  detector, 
where  the  field  is  maximum.  If  the  disk  were  replaced  with  a  focussed  transducer 
the  maximum  coherent  echo  would  be  expected  to  come  from  the  focal  region  which 
is  at  a  distance  from  the  tranducers. 

We  see  also  that  the  decay  time  constant  of  Eg(t),  a/a  *  (5000)"^  secs  is 
just  12  times  longer  than  the  sound  transit  time  2a/c  across  a  disk  transducer 
diameter.  Thus  to  maintain  the  frequency  of  this  signal  high  enough  to  be 
within  the  passband  of  the  transducers  used,  the  diameter  of  the  disk  must  be 
kept  as  small  as  possible.  Since  a  point  source  does  not  have  a  good  low  fre¬ 
quency  response,  this  should  be  used  as  the  transmitter  rather  than  the 
receiver.  Using  a  focussed  element  instead  of  the  disk  will  probably  also  help 
to  increase  the  frequency  of  E5(t). 


III.  Discussion 

As  pointed  out  earlier,  since  neither  an  electro-magnetic  nor  an  acoustic 
transducer  can  transmit  a  d.c.  signal,  the  time  integral  of  all  transmitted 
electromagnetic  or  acoustic  signal  amplitudes  is  zero.  Using  this  fact,  we 
showed  that  the  echoes  of  such  signals  from  constant  density  random  scatterers, 
with  attenuation  neglected,  are  purely  incoherent  for  either  a  plane  wave  or  a 
point  source/emitter/receiver  (i.e.,  E(t)  has  a  time  average  of  zero  when  the 
scatterers  move  with  respect  to  the  transducer).  In  the  presence  of  strong 
field  gradients  of  sufficiently  low  symmetry,  however,  such  as  those  produced 
by  tranducer  pairs,  we  predict  the  existence  of  coherent  echo  components,  de¬ 
fined  as  echoes  with  a  finite  ensemble  average. 

As  shown  by  eqt.  (7)  the  ratio  between  the  coherent  and  incoherent  echo 
powers  is 

Pcoh  .  N  ~g2|  fG(r,t)dvl2 
pinc  ^/|G(r,  t)  | ^dv 

If  the  coherent  echo  comes  mainly  from  a  limited  region  such  as  a  transducer 
focus,  then  a  measurement  of  the  ratio  Pcoh/pinc  should  enable  the  local  value 
of  Np  /p^~i to  be  estimated  when  the  field  function  is  known.  This  should  be 
applicable  to  the  estimation  of  grain  size  in  metals,  and  possibly  to  tissue 
characterization. 


8.  Echo  E^(t)  received  from  a  uniform  scattering  medium  by  the 
point  transducer  pair  of  Figure  7  for  impulse  acceleration  of  the 
surface  at  time  zero. 


Metallurgical  Sample  Preparation 


SPECIMEN  PREPARATION 


Test  specimens  suitable  for  this  investigation  should  be  capable  of  render¬ 
ing  an  unambiguous  ultrasonic  signal  from  fine,  uniformly  distributed  microcon¬ 
stituents.  A  graded  series  of  specimens  should  contain  such  entities  in  sizes 
ranging  in  diameter  from  25-150  jan.*  Physically,  the  specimen  test  faces 
should  be  parallel  and  spaced  not  less  than  50  mm  apart. 

HIGH  PURITY  IRON 

A  hot-forged  bar  of  high  purity  iron  was  selected  as  the  first  candidate 
material  for  test  specimens.  The  bar  was  large  enough  (about  75  mm  In  diameter) 
and,  by  virtue  of  its  composition,  should  have  had  a  microstructure  consisting 
only  of  equiaxed  grains.  That  in  fact  turned  out  to  be  so,  with  one  important 
exception:  the  grain  size  of  the  iron,  shown  in  Figure  9,  was  huge  and  much  too 
large  for  the  Intended  use.  Nevertheless  the  attributes  of  this  material  made 
it  worthwhile  to  see  if  it  could  be  altered  and  saved. 

The  ferrite  grain  size  of  iron  may  be  modified  -  increased  or  decreased  - 
through  appropriate  heat  treatment.  The  heat  treatment  consists  of  repeated 
cycles  of  rapid  heating  to  35°C  above  the  austenitic  transformation  temperature 
(910°C),  holding  until  the  center  of  the  specimen  reaches  temperature,  then 
quenching  in  water  to  reform  ferrite  grains.  This  treatment  sequence  was 
applied  to  one-centimeter  cubes  ten  times,  with  the  results  shown  in  Figure  10. 
Clearly  a  large  reduction  in  grain  size  was  achieved;  however,  the  grains  are 
not  of  uniform  size,  and  some  remain  quite  large.  It  is  questionable  therefore 
whether  specimens  having  such  a  microstructure  would  be  suitable  for  the  purpose 
of  this  investigation. 


*Equiaxed  metallic  grains  in  this  range  of  diameters  correspond  to  standard 
ASTM  grain  size  numbers  between  7.5  and  2.5.  Refer  to  the  attached  Table  of 
Micro-Grain  Size  and  Grain  Size  Charts. 


While  work  on  modifying  the  grain  structure  of  the  high-purity  iron  was  in 
progress,  a  quantity  of  commercially  pure  (i.e.,  1100)  aluminum  was  obtained 
from  Alcoa.  The  ingot  was  cast  by  the  D.C.  process  and  refined  with  a  titanium- 
boron  grain  refining  addition.  Commercial  purity  (1100)  aluminum  contains  small 
quantities  of  iron  and  silicon  as  the  principal  impurity  elements.  Iron/sllicon 
precipitates  appear  as  black  particles  in  the  dendritic  microstructure  shown  in 
Figure  11.  This  material  may  prove  to  be  useful  for  making  specimens;  its  micro¬ 
constituents  are  closely  spaced  and  uniformly  distributed. 

Contingent  Material 

In  the  event  that  neither  the  commercially  pure  aluminum  nor  high-purity 
iron  yields  satisfactory  specimen  material,  it  will  become  necessary  to  place  a 
special  order.  That  situation  has  the  advantage  of  our  being  able  to  specify 
precisely  both  the  physical  dimensions  and  requisite  microstructure  of  the 


material  as  conditions  of  purchase. 


Figure  10.  Pure  Iron  After 
Cyclic  Grain  Refining  Heat 
Treatment.  100X 


Figure  11.  DC  Cast  1100 
Aluminum  Grain  Refined. 
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0.040 

1.58 

0.0355 

1.40 

IX 

248 

15  600 

625 

403 

6.5 

0.038 

1.49 

0.0335 

1.32 

1.43 

221 

18  600 

701 

453 

0.035 

1.38 

0.0311 

1.22 

IX 

IX 

23  000 

816 

32.7 

7.0 

0.032 

1.25 

0.0282 

1.11 

1.01 

IX 

31  OX 

992 

64.0 

0.030 

1.18 

0.0287 

1.05 

OX 

1.40 

37  OX 

1  no 

71.7 

7.3 

0.027 

1.05 

0.0237 

0.933 

0.713 

1.10 

52  500 

1  400 

90.5 

0.025 

0.964 

0.0222 

0.874 

0.825 

0.969 

64  OX 

i  ex 

103 

8.0 

0.0224 

0.884 

0.0199 

0.785 

0X4 

0.781 

88  400 

I  980 

128 

0.0200 

0.787 

0.0178 

0.899 

0.40 

0.6X 

125  OX 

25X 

161 

65 

00189 

0.743 

0.0168 

0.680 

0X6 

0552 

149  OX 

2  810 

181 

9.0 

0.0159 

0.625 

0.0141 

0.555 

0X2 

0X1 

250  OX 

3  970 

256 

0.0150 

0.591 

0.0133 

0.524 

0X5 

0249 

296  OX 

4  440 

287 

9.5 

0.0134 

0.326 

0.0119 

0.467 

0.178 

0276 

420  OX 

5  610 

362 

10.0 

0.0112 

0.442 

0.00997 

0.392 

0.1X 

0.195 

707  OX 

7  940 

512 

0.0100 

0.394 

0.00888 

0.350 

0.10 

0.155 

l.Xx  10“ 

10  OX 

645 

10.5 

000844 

0.372 

0.00838 

0.3X 

0.089 

0.138 

1.19x10“ 

11  2X 

724 

000900 

0.354 

0.00799 

0.315 

0.081 

0.126 

137x10“ 

12  3X 

797 

0.00600 

0.315 

0.00710 

0.280 

0.064 

0.0992 

1.95x10“ 

is  ex 

I  010 

11.0 

000794 

0.313 

0.00705 

0.278 

0.063 

0.0977 

2.Xx  10" 

15  9X 

1  020 

000700 

0.276 

0.00622 

0.245 

0.049 

0.0760 

2.92  x  10“ 

20  4X 

1  320 

11.5 

000667 

0.263 

0.00503 

0.233 

0.045 

0.0691 

336 x  10“ 

22  400 

1  450 

000600 

0.236 

0.00533 

0.210 

0X6 

0.0558 

4.63x  10” 

27  800 

1  790 

12.0 

0.00661 

0.221 

0.00498 

0.196 

0X1 

0.0488 

5.66  x  10“ 

31  7X 

2060 

0.00600 

0.197 

0.00444 

0.175 

0.025 

0.0388 

8.XX  101 

40  OX 

2  580 

12.5 

000472 

0.186 

0.00419 

0.165 

0.022 

0.0345 

951x10“ 

44  900 

2  900 

0.00400 

0.158 

0.00355 

0.140 

0.0160 

0.0248 

15.62X10“ 

62  5X 

4  030 

iao 

0.00097 

0156 

0.00352 

0.1X 

0.0158 

0.0244 

16.0  X  10“ 

63  500 

4  IX 

las 

0.00034 

0.131 

0.00296 

0.117 

0.011 

0.0173 

26.9  x  10* 

89  800 

5  800 

0.00000 

0.118 

0.00266 

0.105 

0.009 

0.0140 

37.0  x  10“ 

111  OX 

7  170 

14.0 

0.00281 

0.111 

0.00249 

0.0981 

0.0079 

0.0122 

45.2  x!0“ 

127  OX 

82X 

0.00250 

0.096 

0.00222 

0.0874 

0.00625 

0.00969 

64.0  X 10“ 

160  OX 

10  3X 

'The  use  of  X  is  recommended  instead  of " — 1”  or  "minus  1"  to  avoid  confusion. 


'Value  of  Heyn  intercept  for  equfaxed  grains. 

From:  Designation  E-112,  A  STM  Standards,  1961,  Part  3—  Metal  Test  Methods;  published  by  American  Society  for  Testing  and 
Materials^hiladelphia,  Pa. 
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then  e(t)  can  be  written  as 
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where  e'(t)  is  Independent  of  a 

3,2 

When  t'  ♦o,  e'(t')  ♦  £_£ _ 
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3  » 2 

here  c  is  the  echo  from  scatterers  uniformly  distributed  in  a  hemisphere  due 
4 

to  an  Impulse  exerted  on  a  point-transmitter/infinite  plane-reviewer  transducer 
pair 

when  t  ’  «  ,  e’(t')  >  £_ 

4  IT 


here  £_  is  the  echo  from  scatterers  uniformly  distributed  in  half  space  due 
4  n 

to  an  impulse  exerted  on  a  point-transmitter/point-receiver  transducer  pair. 
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Several  computer  programs  are  attached* 


EdGltal .  f 


This  program  is  used  for  calculating  the  impulse 
response  of  a  disc-point  transducer  pair 
for  the  case  when  the  time  t  <  a/c 

where  a  is  the  radius  of  the  disc  transducer  and  c  the 
velocity  of  sound  in  the  medium  .  The  program  uses 
subroutine  DBLIN  of  the  IMSL  library. 


INTEGER  IER 

real  DBLIN, Fl, AX, AY,BX, BY, AERR, ERROR, Cl 

EXTERNAL  Fl , AY, BY 

common  /CI/a,c,t 

real  ar ( 400 ) , swt ( 400 ) 

character  *10  file 

Parameter  a  is  the  radius  of  the  disc  transducer,  c  is  the  veloci 
of  sound  in  the  medium  and  AERR  the  absolute  error  of  the 
calculation  by  the  definition  of  the  DBLIN  subroutine. 

The  final  calculated  value  of  the  impulse  response  is  stored  in 
in  an  array  named  as  ar  while  discrete  tine  is  stored  in  array  sw 


When  compiling  &  linking,  type 

ill  -o  czc  Edeltal.f  /usr/lib/imsla. a 


a-1 . 25 
c-150000. 0 
AERR-0 .0001 

write (6 , 10 ) 

format (4x, ' Input  t  as  starting  time:'} 
read(5,20)  t 
format (f 16. 12) 
write (6,40) 

format ( 2x, ' Input  the  name  of  the  file  to  store  the  calculated  dat 

read (5.60)  file 

format  f alO ) 

t=0. 000004893047 

open ( unit-1 , file- ' dataEE22 1 , status= ' new ' , access- ' sequential ' , 

form- ' formatted 1 ) 

do  150  km=l ,10 

t-t*l . 05 

tnr-a/c 

if(t  .gt.  thr  .or.  t  .eq.  thr)  go  to  111 

AX-0.0 

BX  -a 

Cl -DBLIN ( Fl , AX, BX, AY, BY, AERR, ERROR, IER) 
sw  t  ( km )  - 1 
ar ( km) =C1 

wr i te ( 1 , 2  02 ) swt ( km) , ar ( km) 


format (f 16 . 12 , fl6 . 8 ) 

write(6,126)  km, swt ( km) , ar (kn) 

format (lOx, i5 ,3x, fl6. 12, 3x, £16. 8 ) 

continue 

close ( unit=l ) 

stop 

END 

REAL  FUNCTION  Fl(x,y) 
real  x,y 
common  /CI/a,c,t 

To  avoid  the  difficulty  due  to  the  inability  of  computer  to  deal 
with  the  case  where  0/0  appears  in  the  calculation,  we  ignore  a  small 
part  of  scatterers  near  the  origin. 

rr=sqrt  (x*x-;-y*y ) 

R-sqrt  (y*y-:-  (a-:-x)  *  (a-:-x)  ) 

Rp-sqrt  (y*y-:-  (a-x)  *  (a-x) ) 

u-1 0000000000 .0* (c*t-rr) * (c*t-rr ) 

u-- u-1 00 00000000 . 0*y*y-rl 000 000 00 00 . 0*x*x-10000000000 . 0*a*a 

v-2. 0*x*100000.0*sqrt(10000000000.0*(c*t-rr)*(c*t-rr)-l 0000000000. 0*y*y) 

if(rr  .It.  0.000125)go  to  598 

ind-0 

j  nd-0 

iffrr  .gt.  (c*t-R)  .and.  rr  .It.  (c*t-Rp))  ind-5 

i f ( i nd  .eq.  5)  Fl-100000 .0*x*c/ (2.0*3 . 1416*3 . 1416*100000 . 0*rr) *acos ( u/v) 

iffrr  .gt.  (c*t-y)  .or.  rr  .eq.  (c*t-y) )  Fl-0.0 

if(rr  .gt.  (c*t-Rp)  .and.  rr  .It.  (c*t-y) )  jnd-3 

iffjnd  .eq.  3)  Fl=100000 . 0*x*c/ (2.0*3. 1416*rr*100000.0) 

iffrr  .eq.  (c*t-Rp) )  Fl=100000.0*x*c/(2.0*3.1416*rr*100000.0) 

iffrr  .It.  (c*t-R)  .or.  rr  .eq.  (c*t-R) )  Fl-0.0 

go  to  577 

Fl-0 . 0 

RETURN 

END 


REAL  FUNCTION  AY(x) 

real  x 

AY=0 . 0 

return 

end 

REAL  FUNCTION  3Y(x) 
real  x 

common  /CI/a,c,t 

3Y-a/2 . 0 
return 


Edel ta2 . f 


This  program  is  used  for  calculating  the  impulse 
response  of  a  disk-point  transducer  pair 
for  the  case  v/hen  the  time  t  >  a/c 

where  a  is  the  radius  of  the  disc  transducer  and  c  the 
velocity  of  sound  in  the  medium  .  The  program  uses 
subroutine  DBLIN  of  the  IHSL  library. 

INTEGER  IER 

doubl e  precision  DBLIN, FI , F2 , AX, AY, BX, BY, AAY, BBY, AERR, ERROR, Cl , C2 
EXTERNAL  Fl , F2 , AY, BY, AAY, BBY 
common  /CI/a,c,t 

double  precision  ar (400) ,swt{400) 
character  *12  file 

Parameter  a  is  the  radius  of  the  disc  transducer,  c  is  the  velocity 
of  sound  in  the  medium  and  AERR  tne  absolute  error  of  the 
calculation  by  the  definition  of  the  DBLIN  subroutine.  Since  the 
calculated  value  of  impulse  response  for  t>  a/c  is  larger, 
when  setting  AERR  equal  to  10.0,  the  relative  error  is  still 
small  enough  to  be  accepted.  The  final  calculated  value  of 
the  impulse  reponse  is  stored  in  an  array  named  as  ar  while  discret 
time  is  stored  in  array  swt. 


Nhen  compiling  &  linking,  type 

r'77  -o  zcz  Edelta2.f  /usr/lib/imsls.  a 

a-0 . 125 
c=150000.0 
AERR=0 . 1 


write(6,10) 

format  ( 4>:, '  Input  t  as  starting  time:') 
read  (5,20)  t 
format (fl6.12) 
write(6 ,40) 

format (2x, ' Input  the  name  of  the  file  to  store  the  calculated  daea: 

read(5,50)  file 

format (al2) 

t==0. 0000701152661 

open ( unit-1 ,  f  ile=  ' datasal ' ,  status- '  new'  ,  access--  *  sequential '  , 
f orm= 1  formatted' ) 
do  150  km=l , 20 
t==t*1.02 

a:l=o.o 

BX-a 

Cl-DBLIN (Fl , AX, BX, AY, BY, AERR, ERROR, IER) 


BX~c*t/2 .  O-i-a/2 . 0 
sw  t ( km ) = t 

C2=DBLIN(F2,AX,BX, AAY,BBY, AERR, ERROR, IER) 

ar (kra)=Cl+C2 

AERR=abs (ar ( km) *0.01) 

write (1, 202 )swt (km) , ar (km) 

format (fl6.12,fl6.8) 

write(6,126)  km,swt(km) ,ar (km) 

format (10x,i5,3x,fl6.12,3x,fl6.8) 

continue 

close ( unit=l ) 

stop 

end" 


FUNCTION  Fl(x,y) 
double  precision  x,v,Fl 
common  /CI/a,c,t 

To  avoid  the  difficulty  due  to  the  inability  of  computer  to  deal 
with  the  case  where  0/0  appears  in  the  calculation,  we  ignore  a  small 
part  of  scatterers  near  the  origin. 

r  r  =  sq r  t  ( x *x-.-y  *y ) 

R-sqrt (y*y+(a+x) *  (a+x) ) 

Rp~sqrt  (y*yJr(a-x)  *(a-x) ) 

u-1 000 00000 . 0* (c*t-rr) *(c*t-rr) 

u-u-100000000 .  Q*y  *y-;-100000000 . 0*x*x-100000000 .0*a*a 

v=2.0*x*10000.0*sqrt (100000000. 0* (c*t-rr ) * (c*t-r r) -100000000. 0*y*v) 

if ( rr  .It.  0 .0000125 ) go  to  598 

ind-0 

jnd-0 

if(rr  .gt.  (c*t-R)  .and.  rr  .It.  (c*t-Rp))  ind-5 

if ( ind  .eq.  5)  Fl=10000. 0*x*c/(2. 0*3. 1416*3. 1416*10000. 0*rr) *acos(u/v) 

if(rr  .gt.  (c*t-y)  .or.  rr  .eq.  (c*t-y))  Fl-0.0 

if(rr  .gt.  (c*t-Rp)  .and.  rr  .It.  (c*t-y))  jnd=3 

if ( jnd  .eq.  3)  Fl=10000.0*x*c/(2.0*3.1416*rr*10000.0) 

if ( rr  .eq.  (c*t-Rp))  F1^10000 . 0*x*c/ (2 . 0*3 . 1416*rr*10000 . 0 ) 

if(rr  .It.  (c*t-R)  .or.  rr  .eq.  (c*t-R))  Fl=0.0 

go  to  577 

F1=0 . 0 

RETURN 

END 


FUNCTION  F2(x,y) 
double  precision  x,y,F2 
common  /CI/a,c,t 

rr=sqrt(x*x-:-y*y) 

R-sqrt (y*y+ (a+x) * (a+x) ) 

Rp-sart(y*y-i-(a-x)  *(a-x) ) 

u=l 000 00 00 0.0* (c*t-rr) * (c*t-rr) 

u«  u-100000000. 0*y*y:-l  00000000. 0*x*x-100000000.0*a*a 


v=2.0*x*10000.0*sqrt( 100000000. 0* (c*t-rr ) * (c*t-rr) -100000000. 0*y*y) 
knd=0 

if(rr  .gt.  c*t-R  .and.  rr  .It.  c*t-Rp)  knd=8 

if ( knd  .eq.  8)  F2^10000.0*x*c/<2. 0*3. 1416*3. 1416*rr*10000.0)*acos(u/v) 

if(rr  .gt.  (c*t-Rp)  .or.  rr  .eq.  (c*t-Rp))  F2-0.0 

if ( r r  .It.  { c*t-R)  .or.  rr  .eq.  (c*t-R))  F2-0.0 

RETURN 

END 

FUNCTION  AY ( X ) 

double  precision  x,AY 

AY-0 . 0 

return 

end 

FUNCTION  R Y ( X ) 
double  precision  x,BY 
common  /CI/a,c,t 

BY-c*t/2 . 0 

return 

end 


FUNCTION  AAY(x) 
double  precision  x, AAY 

AAY-0 . 0 

return 

end 

FUNCTION  BBY(x) 
double  precision  x, BEY 
common  /CI/a,c,t 

BBY-c*t/2.0* (1.0-a*a/ (c*c*t*t) ) 
return 
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1030 


1045 

1055 


1205 


joindata.  f 


Programs  Edelta2  or  Edoltal  might  stop  being  executed  while 
they  are  supposed  to  be.  The  reason  can  be  the  singularity 
or  ill  behaviour (not  smooth  enough  according  to  the  request  of 
given  error)  of  the  integrand.  The  easy  way  to  get  rid  of 
such  difficulty  is  to  jump  over  these  points  where  the 
integrand  doesn't  behave  well.  It  will  not  do  harm  to  our 
calculation  because  such  points  are  not  continuous.  From  the 
point  of  view  of  physics,  the  calculated  curve  should  be 
continuous.  This  Program  joindata. f  is  used  to  join  separated 
data  files  together  to  produce  a  final  data  file. 

real  swt (1000) ,ar (1000) 
character*8  file 
character*!  q,aq,qqq,q4 ,q5 


nor=l 

write(6,1000) 

f ormat (2x, 1  Do  you  want  to  read  some  data  file  ?  y/n?') 
read(5,1010)  q 
format (al) 

if (q  .eq.  'n')  go  to  1045 
write(6,ll) 

format (2 x, ' Input  the  name  of  the  file:') 
read(5,14)  file 
format (a8 ) 
write(6, 1020) 

f ormat (2x, ' How  many  records  does  the  file  have?') 
read(5,1030)  norl 
nom=nor+norl-l 
format ( i4 ) 

open ( unit=l , file- file, status= 'old' ,access= ' sequential 1 , 

f orm= 'formatted' ) 

do  100  i=nor,nom 

read ( 1,15) swt (i) ,ar<i) 

format (f 16.8, f 16. 8) 

write(6,10) i,  swt(i),ar(i) 

format(4x,  i4 , f 16 . 8 , f 16 . 8 ) 

continue 

close (unit=l) 

nor=nom+l 

go  to  77 

write(6,1055) 

format (3x, ’ Do  you  want  to  type  in  some  data?  y/n?') 

read(5,1205)q5 

format (al ) 

if(q5  .eq.  'n')  go  to  1095 
nor-nor-1 


write (6,1200) 

format (3x, ' will  t  automatically  be  calculated?  y/n?') 
read(5,1210)  qqq 
f ornat (al ) 

if (qqq  .eq.  'n')  go  to  1250 
write (6,1212) 

format(2x, ' Is  this  the  first  value  of  t?  v/n?') 
read(5,1215)  q4 
format (al ) 

if(q4  .eq.  'n')  go  to  1265 
write(6,1260) 

format (2x, ' Input  previous  t  as  initial  value:1) 
read(5,1270)  t 

if (q4  .eq.  'y')  t=t*l . 002*1 . 002 
nor=nor+l 
format (f 16 . 12 ) 
t=t*l. 002 
write(6f1275) t 
format(fl6.8) 
write(6,1280) 

format (4x, ' Input  echo(t):') 
read(5,1290)  ar(nor) 
format (fl6. 8) 
swt (nor) =t 
go  to  108 
nor=nor-:*l 
write(6,1060) 

format (2x, ' Input  t  &  echo(t):') 
read(5,1070)  swt (nor ) ,ar (nor) 
f ormat (f 16 . 12 , f 16 . 8 ) 
write(6,1080) 

format(2x, 'Do  you  think  you  already  input  all  data?  y/n?') 
read  (5,1090)  qq 
format (al ) 

if (qq  .eq.  'n')  go  to  1050 
nor-nor-1 
write(6, 1100) nor 

format (5x, 'The  Mo.  of  records  of  the  new  file  is  :',i5) 
open  ( unit=l ,  f  ile=  '  CTFl2data ' ,  status=  '  new'  , access-”  '  sequential ' 
1  form= ' formatted' ) 
do  1150  j=l,nor 
write(l,1110)swt(j) ,ar(j) 
format (f 16.8, fl6. 8) 
continue 
op 


planepoint . f 


This  program  is  used  to  calculate  the  echo  from  scatterers 
which  are  uniformly  distributed  in  half  space  due  to  the 
impulse  input  to  infinite  plane-point  transducer  pair.  The 
calculated  data  can  serve  as  a  comparison  with  the  calculated 
result  by  program  Edeltal.fi.  Note  the  data  is  normalized  here 
by  being  divided  by  2.0*3.1416. 

character*6  file 


write(6,ll) 

format (2x, ' input  the  name  of  the  file:') 
read(5 , 14) file 
format (a6 ) 

open  (uni t=l ,  f  ile=f  ile,  status=  '  new'  ,  access^ '  sequential ' , 

f orm= ' formatted' ) 

t^O. 0000001 

do  100  i=l ,450 

t=1.01*t 

f =150000. 0 *150000. 0*150000. 0*t*t/ (4. 0*2. 0*3. 1416) 

write(l,20)  t,f 

format(fl6.12,fl6.8) 

continue 

close(unit=l) 


curvef it. f 


real  swt (1000) ,ar(1000),aa(1000),rna(5) 
character*12  f ile, f ilel ,f ile2 


write(6,30) 

format (4x, ' Input  the  name  of  the  file:1) 
read(5,40)  file 
format (al2 ) 

open(unit=l,f  ile=f  ile,  status-- '  old*  , access= ' sequential '  , 

forn= ' formatted' ) 

do  10  i=l,387 

read  (1,20)  swt(i),ar(i) 

format (£16 . 12, fl6 . 8 ) 

continue 

close(unit=l) 

write (6 ,70) 

f ormat (2x, ' Input  the  name  of  the  new  file  to  store  the  fitting  curve:') 
read(5,80)  filel 
format (al 2) 

write(6,72) 

format (4x, ' Input  the  name  of  the  new  file  to  store  the  estimates:') 
read(5,80)  £ile2 

Input  the  initial  value  of  amp,  sigma  &  cst  as  curve  fitting  parameters 
write(6,50) 

format (4x, ' input  the  initial  value  for  amp, sigma  6  cst  :') 
read(5,60)  amp, sigma, cst 
f ormat ( f 20. 10, f 20. 10, f 20. 10) 

de=l . 0 
all=l .  0 
rin-0 .001 


curve  fitting  for  estimating  amp, sigma  &  cst 

call  difference (amp, sigma,cst,swt,ar,e,aa) 
ampl= amp-0 .01 
sigmal=sigma-0 .01 
cstl-cst-0.01 

call  dif f erence (arapl, sigma, cst,swt,ar,ea,aa) 
call  dif f erence (amp, sigmal , cst, swt , ar , es, aa) 
call  difference (amp, sigma, cstl , swt,ar, ec, aa) 
ga= (e-ea)/0.01 
gs= (e-es)/0.01 
gc= ( e-ec) /0 . 01 


adl®amp-all*ga 

sdl=sigma-all*gs 

cdl=cst-all*gc 

call  dif ference(adl, sdl,cdl, swt, ar , el,aa) 

if  (e  .It.  el)  go  to  390 

al2-all*1.3 

ad2=anp-al2*ga 

sd2=sigma-al2*gs 

cd2=cst-al2*gc 

call  dif f erence (ad2, sd2 , cd2 , swt , ar, e2 , aa) 

if(e2  .gt.  el)  go  to  360 

el=e2 

all=al2 

al2-all*1.3 

ad2=anp-al2*ga 

sd2=sigma-al2*gs 

cd2=cst-al2*gc 

call  dif f erence (ad2 , sd2 , cd2, swt, ar ,  e2  ,  aa) 
go  to  370 

e2=el 

al2-all 

all®0.618*al2 

adl=amp-all*ga 

sdl=sigma-all*gs 

cdl=cst-all*gc 

call  difference (adl, sdl, cdl, swt,arf  el ,  aa) 
iffel  .gt.  e)  go  to  390 

de-abs ( e-el ) 
rna(l) =rna(2) 
rna(2)=rnaf3) 
rna(3)=rna(4) 
rna(5)=de 

srna-rna(l)  -;-rna(2)  +rna(3)  +rna  (4)+rna  (5) 

if(srna  .It.  rin)  go  to  130 

e=el 

amp=adl 

sigma- sdl 

cst-cdl 

go  to  380 

e2=el 

al2=all 

all=al2*0.618 

adl=amp-all*ga 

sdl®sigraa-all*gs 

cdl=cst-all*gc 

call  dif f erence (adl , sdl , cdl, swt , ar , el , aa) 
go  to  400 

open ( unit=l , f ile=f ile2 , status® 1  new' , access® 'sequential ' 
1  form® ' formatted' ) 


write (1,160)  amp, sigma, cst 
format ( f 20 .12) 
close( unit-1 ) 

Create  new  data  file  to  be  displied 

open ( unit=l , f ile=f ilel , status= ' new ' , access= ' sequential ' , 

form- 1  formatted' ) 

do  90  jmm-~l,387 

write(l,100) swt( jmm) ,aa(jmm) 

format (fl6.12,fl6.8) 

continue 

close ( unit=l ) 

write ( 6 , 137 ) 
format (2 x, ' The  end. ' ) 

stop 

end 


subroutine  difference  (AI1P,  SIGNA,  CST,  SUT,  AR,  E,  AA) 
REAL  AR(IOOO) ,AA( 1000) ,SNT( 1000) 

E=0 . 0 

DO  200  KK=1 , 3 87 

AA  ( KK )  =  AI1P*EXP  ( -1 . 0  *S IGriA*SITT  ( KK)  )  +CST 
E=  E+ABS  ( AA(KK)  -AR(  KK)  ) 

CONTINUE 

RETURN 

END 
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